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Abstract. Given a birth-death process on N with semigroup (Pt)t>o and a discrete 
gradient d u depending on a positive weight u, we establish intertwining relations of 
the form d u Pt = Qtd u , where {Qt)t>o is the Feynman-Kac semigroup with potential 
V u of another birth-death process. We provide applications when V u is non-negative 
and uniformly bounded from below, including Lipschitz contraction and Wasserstein 
curvature, various functional inequalities, and stochastic orderings. Our analysis is 
naturally connected to the previous works of Caputo-Dai Pra-Posta and of Chen on 
birth-death processes. The proofs are remarkably simple and rely on interpolation, 
commutation, and convexity. 



1. Introduction 

Commutation relations and convexity are useful tools for the fine analysis of Markov 
diffusion semigroups [B-E\ IB"} [L] . The situation is more delicate on discrete spaces, due 
to the lack of a chain rule formula (EU \M ICheTl U-P\ UTTl ICha2l IC-DP-P1 IChe3] . In 
this work, we obtain new intertwining and sub-commutation relations for a class of birth- 
death processes involving a discrete gradient and an auxiliary Feynman-Kac semigroup. 
We also provide various applications of these relations. Our analysis is naturally related 
to the curvature condition of Caputo-Dai Pra-Posta |C-DP-P] and to the Chen exponent 
of Chen |CheltlChe3] . More precisely, let us consider a birth-death process (X t ) t > on the 



state space N := {0, 1,2,.. .}, i.e. a Markov process with transition probabilities given by 

(\ x t + o(t) ify = x + l, 

P t x (y) = F x (X t = y) = L x t + o(t) if y = x - 1, 

|l - (Aa, + u x )t + o(i) if y = x, 

where \im t ^ot~ 1 o(t) — 0. The transition rates A and v are respectively called the birth 
and death rates of the process (X t )t>o- The process is irreducible, positive recurrent (or 
ergodic), and non-explosive when the rates satisfy to A > on N and v > on N* and 
vq = and 



< oo and ^ — + — + • • • + -— 

x=l v \ v 2''' v x x=1 \A X A x A x _i Ax'''M*C 



OO, 



respectively. In this case the unique stationary distribution \x of the process is reversible 
and is given by 

^)=/i(0)II— ,XGM with /x(0):= l + . (1.1) 

y=l V V V x=l ^1^2 -~Vx j 



Date: Preprint accepted for publication in Bernoulli. 
2000 Mathematics Subject Classification. 60J80 ; 39C05. 

Key words and phrases. Birth-death process; Feynman-Kac semigroup; discrete gradients; intertwining 
relation; functional inequalities. 

1 



2 



DJALIL CHAFAI AND ALDERIC JOULIN 



Let us denote by T (respectively T+ and Tj) the space of real- valued (respectively positive 
and non-negative non-decreasing) functions / on N. The associated semigroup {Pt)t>o is 
defined for any bounded or non-negative function / as 

oo 

p t f(x) = E x [f(x t )) = x: f(y)P t x (y), x eN. 

This family of operators is positivity preserving and contractive on L p (p), p G [l,oo]. 
Moreover, the semigroup is also symmetric in £ 2 (p) since X x fi(x) = ui +x fj,(l + x) for any 
iGN (detailed balance equation). The generator C of the process is given for any / G J 
and iGNby 

Cf(x) = X x (f(x + 1) - f(x)) + v x (f(x - 1) - f(x)) 

= X x df(x) + u x d*f(x), 

where d and d* are respectively the forward and backward discrete gradients on N: 

df(x):=f(x + l)-f(x) and d* f(x) := f(x - 1) - f(x). 

Our approach is inspired from the remarkable properties of two special birth-death pro- 
cesses: the M/M/l and the M/M/oo queues. The M/M/oo queue has rates \ x = A 
and u x = vx for positive constants A and v. It is positive recurrent and its stationary 
distribution is the Poisson measure p p with mean p = A/p. If 3$ x # stands for the bino- 
mial distribution of size x G N and parameter p G [0, 1], the M/M/oo process satisfies 
for every igN and t > to the Mehler type formula 

Sf(X t \X = x) = SS x ^t * Hp^_ e -vty (1.2) 

The M/M/l queening process has rates A^ = A and v x = i/l N \{ } where < A < v 
are constants. It is a positive recurrent random walk on N reflected at 0. Its stationary 
distribution p is the geometric measure with parameter p := \ jv given by p(x) = (1— p)p x 
for all x G N. A remarkable common property shared by the M/M/l and M/M/oo 
processes is the intertwining relation 

dC = C v d (1.3) 
where C v = £ — V is the discrete Schrodinger operator with potential V given by 

• V(x) := v in the case of the M/M/oo queue 

• V(x) := z/l {0} (x) for the M/M/l queue. 

Since V > in these two cases, the operator C v is the generator of a birth-death process 
with killing rate V and the associated Feynman-Kac semigroup (P^) t >o is given by 



pVf( x ) = E x 



/(X t )exp(- J*V(X s )ds 



The intertwining relation (II .3|) is the infinitesimal version at time t = of the semigroup 
intertwining 



dP t f(x) = P t v df(x)=E. 



. t 

V; 



df(X t ) exp (-J V(X s )ds 



;i.4) 



Conversely, one may deduce (11.41) from fll.3p by using a semigroup interpolation. Namely, 
if we consider s G [0,t] !->■ J(s) := PjOPt^sf with V as above, then (11.41) rewrites as 
J(0) = J(t) and ([I3D follows from (TOD since 

J'(s) = Pj (£ v dP t _J - dCPtsf) = 0. 



In section [2], we obtain by using semigroup interpolation an intertwining relation similar 
to fll.4p for more general birth-death processes. By using convexity as an additional 
ingredient, we also obtain sub-commutation relations. These results are new and have 



INTERTWINING AND COMMUTATION RELATIONS FOR BIRTH-DEATH PROCESSES 3 



several applications explored in section [3J including Lipschitz contraction and Wasserstein 
curvature (section |3.1|) . functional inequalities including Poincare, entropic, isoperimetric 
and transportation-information inequalities (section [3.2|) . hitting time of the origin for the 
M/M/l queue (section 13.31) . convex domination and stochastic orderings (section 13.41) . 

2. Intertwining relations and sub-commutations 

Let us fix some u G J±. The w-modification of the original process (X t )t>o is a birth- 
death process (X U)t )t>o with semigroup (P u ,t)t>o and generator C u given by 

£ u f(x) = \ u x df(x) + v u x d*f(x), 

where the birth and death rates are respectively given by 

K '■= K+i and vl := v x . 

One can check that the measure Aw 2 /i is symmetric for (X u<t )t>o- As consequence, the 
process (X U)t )t>o is positive recurrent if and only if Aw 2 is /z-integrable. From now on, we 
restrict to the minimal solution corresponding to the forward and backward Kolmogorov 
equations given as follows: for any function / G T with finite support and t > 0, 

r, Pu,tf Pu,t£"uf ^uPu.tf) 

at 

cf. |Che2| th. 2.21]. In order to justify in all circumstances the computations present in 
these notes, we need to extend these identities to bounded functions /. Although it is not 
restrictive for the backward equation, the forward equation is more subtle and requires an 
additional integrability assumption. From now on, we always assume that the transition 
rates A" and v u and also the potential V u are P u<t integrable. 
We define the discrete gradient d u and the potential V u by 

d u := (l/u)d and V u {x) := u x+1 - v% + X x - A". 

Let if : R — > R be a smooth convex function such that for some constant c > 0, and for 
all r G R, 

(p'(r)r > cip{r). (2.1) 
In particular, the behavior at infinity is at least polynomial of degree c. 

Let us state our first main result about intertwining and sub- commutation relations 
between the original process (X t ) t > and its u-modification (X u t ) t > . To the knowledge 
of the authors, this result was not known. A connection to Chen's results on birth-death 
processes |Che2] is given in section [3] in the sequel. 

Theorem 2.1 (Intertwining and sub-commutation). Assume that the process is irre- 
ducible, non-explosive and that the potential V u is lower bounded. Let f G T be such that 
su P»/eN \duf{y)\ < oo, and let x G N and t > 0. Then the following intertwining relation 
holds: 

d u P t f(x) = P%d u f(x) = E. 



duf(X Utt ) exp (~ J V u (X UjS )ds 



(2.2) 



Moreover, ifV u >0 then we have the sub-commutation relation 



ip{d u P t f) (x)<E x 



<f(d u f){X u>t ) exp (-J*cV u (X U)S )d s y^ . (2.3) 



Proof. The key point is the following intertwining relation 

d u C = CZ"d u , (2.4) 

where C u is the generator of the ^-modification process (X u j)t>o and CZ« '■= C u — V u 
is the discrete Schrodinger operator with potential V u . Note that the relation (12.41) is 
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somewhat similar to (11. 3p and follows by simple computations. To prove (12.21) from (12.41) . 
we proceed as we did to obtain (II .4|) from (II. 3|) . If we define 

se[o,t]^J{s) := P^; s d u Ptsf, 

then (12. 2p rewrites as J(0) = J(t). Hence it suffices to show that J is constant. By 
[Chelj we know that if d u f is bounded then d u P t - s f is also bounded. Hence using the 
Kolmogorov equations and (12.41) . we obtain 

/(*) = P£ (c^d u P t _J - d u CPt- s f) = 0, 
yielding to the intertwining relation (12. 2|) . 

Now let us prove the sub-commutation relation (12.31) by adapting the previous interpola- 
tion method, under the additional assumption V u > 0. Denoting 

s e [0,t] ^ J c (s) := P%Md u P t - a f), 

then (12. 3p rewrites as J c (0) < J c {t). Hence let us show that J c is a non-decreasing 
function. Since (p(d u P t - s f) is bounded, we have by the Kolmogorov equations: 

J'c(s) = P^(T) where T — C^ u <p(d u P t ^ s f) — <p'{d u P t - a f) d u CP t - s f- 

Letting g u = d u P t - s f, we obtain, by using (12 .4p . 

T = C?M9u) ~ v'(9u)C^g u 

= X u (dip(g u ) - ip'(g u )dg u ) + v u (d*(p(g u ) - (p'(g u )d*g u ) + V u (<p'(g u )g u - op(g u )) 

= X u A^(g u , dg u ) + u u Af{g u , d*g u ) + V u (<p'(g u )g u - op(g u )) 

where A v {r, s) = (p(r + s) — (p(r) — (p'(r)s is the so-called A-transform of (p studied in 
[Cha2] also known in convex analysis as the Bregman divergence associated to ip [BrJ. 
Note that g u + dg u = g u (- + 1) and g u + d*g u = g u (- — 1). Now, since <p is convex, we have 
A v > 0. Moreover, using (12.11) and V u > we obtain that T > 0. Finally, we get the 
desired result since the Feynman-Kac semigroup (P^t u )t>o is positivity preserving. □ 

Remark 2.2 (Ergodic condition). The potential V u in theorem \2.1\ is assumed to be lower 
bounded. When it is positive, the so-called Chen exponent inf^N V u (y) is related to the 
exponential ergodicity of the original process (X t ) t > , cf. |Chel] . However identity (12. 2p 
does not require such an ergodic assumption. A nice study of the exponential decay of 
birth-death processes was recently studied by Chen in [Che3], with special emphasis on 
non- ergodic situations including transient cases. 

Remark 2.3 (Case of equality). According to the proof of theorem \2. 1\ the assumption 
V u > can be dropped if the convex function tp realizes the equality in (12.11) . Such an 
observation was expected since in this case the use of Holder's inequality in (12. 2p entails 
the desired result. 

Remark 2.4 (Propagation of monotonicity) . The identity (12.21) provides a new proof of 
the propagation of monotonicity [HI prop. 4.2.10]: if f G then P t f G for all t > 0. 
See section \3.4\ for an interpretation in terms of stochastic ordering. 

Remark 2.5 (Other gradients). Theorem \2. l\ possesses a natural analogue for the discrete 
backward gradient d* . We ignore if there exists a useful "balanced" intertwining relation 
involving a combination of both forward and backward gradients. 

Remark 2.6 (Higher dimensional spaces). The extension of theorem \2.1\ to higher dimen- 
sional discrete processes such as queuing networks or interacting particles systems arising 
in statistical mechanics is a very natural question, but seems to be technically difficult. 
However a first step has been emphasized by Wu in his study of functional inequalities for 
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Gibbs measures through the Dobrushin uniqueness condition: see step 1 in the proof of 
proposition 2. 5 in [W] . 

Our second new result below complements the previous one for the case u — 1. Let 
X be an open interval of R and let ip : X — > R be a smooth convex function such that 
if" > and —1/ip" is convex on X. Following the notations of |Cha2j . we define on the 
convex subset Ax '■= {(r, s) G R 2 : (r, r + s) G X x X} the non-negative function on 
Ax by 

S v (r, s) := (y'(r + s) — </?'(?")) s, (r, s) G Ax- 

By theorem 4.4 in |Cha2j . B^ is convex on Ax- Some interesting examples of such 
functionals will be given in section 13.21 below. 

Theorem 2.7 (Sub-commutation for 1-modification) . Assume that the process is ir- 
reducible and non- explosive. If the transition rate X is non-increasing and v is non- 
decreasing then for any function f G T such that sup^ gN \df(y)\ < oo and for any t > 0, 

B* (PJ, dPJ) < P$B*(f, df) (2.5) 
where the non-negative potential is V\ := d{y — A). 

Proof. Under our assumption, the two processes (X t ) t > and (Xi )t ) t >o are non-explosive. 
By using standard approximation procedures, one may assume that / has finite support. 
If we define s G [0,t] ^ J(s) := P^\B^(P t _ s f \dP t - 8 f) we see that rewrites as 

J(0) < J(t). Denote F = P t - S f and G = dP t -,f = OF. Since B^(F,G) is bounded, the 
Kolmogorov equations are available and using (12.41) with the constant function u — 1, we 
have J'(s) = P$(T) with 

T = CY 1 B' P (F, G) - -^-B V (F, G) CF - -^-B V (F, G) HY X G 
ox ay 

= X l dB v (F, G) - X^-B^fF, G) OF - X 1 -^-B tp (F, G) dG 
ox oy 

+ u 1 d*B v (F, G) - v-H-B^tF, G) d*F - u l ^-B v (F, G) d*G 
ox oy 

+ d(v - A) {^B^{F, G)G - B V (F, G)j 
> dv (J^B V (F, G)G - B V (F 1 GO J 

- d\ (JLb^F, G)G- ^B v (F, G)G- B*{F, G) \ , 

and where in the last line we used the convexity of the bivariate function B 9 . Moreover, 
since the birth and death rates A and v are respectively non-increasing and non-decreasing 
on the one hand, and using once again convexity on the other hand, we get 

dy v ' ' ~ 1 B*(F,G) 

from which we deduce that T is non-negative and thus J is non-decreasing. □ 

Remark 2.8 (Diffusion case). Actually, the intertwining relations above have their coun- 
terpart in continuous state space, as suggested by the so-called Witten Laplacian method 
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used for the analysis of Langevin-type diffusion processes, see for instance Helffer's book 
[Hj . Let A be the generator of a one-dimensional real-valued diffusion (X t ) t >o of the type 

Af = a 2 f" + bf, 

where f and the two functions a, b are sufficiently smooth. Given a smooth positive 
function a on M, the gradient of interest is V / = a f . Denote (Pt)t>o the associated 
diffusion semigroup. Then it is not hard to adapt to the continuous case the argument of 
theorem \2.1\ to show that the following intertwining relation holds: 

rt 



V a f(X a , t ) exp(-J o V a (X a , s )ds^ 



V a Ptf(x) = E x 

Here (X O)t ) t > is a new diffusion process with generator 

A a f = a 2 f" + b a f 
and drift b a and potential V a given by 



a" . , a' 



b a := 2aa' + b — 2a — and V a :=a b' H b a . 

a a a 

In particular, if the weight a = a, where a is assumed to be positive, then the two processes 
above have the same distribution and by Jensen's inequality, we obtain 



\V a Ptf(x)\<E x 




a 



V a f(X t ) I exp - / aa" -b' + b-\ (X s ) ds 



a 



Hence under the assumption that there exists a constant p such that 



a 



inf aa" — b' + b — > p, 
a 

then we get \V a Ptf\ < e~ pt P t \ V CT /| . This type of sub- commutation relation is at the heart 
of the Bakry-Emery calculus |B-E[ |B| [L] . See also [M-T] for a nice study of functional 
inequalities for the invariant measure under the condition p = 0. However, as we will see 
in remark U75\ below, such a choice of the weight is not really adapted when studying the 
optimal constant in the Poincare inequality. 

3. Applications 
This section is devoted to applications of theorems 12.11 and 12.71 

3.1. Lipschitz contraction and Wasserstein curvature. Theorem 12.11 allows to re- 
cover a result of Chen |Chel] on the contraction property of the semigroup on the space 
of Lipschitz functions. Indeed, the intertwining f!2.2p can be used to derive bounds on the 
Wasserstein curvature of the birth-death process, without using the coupling technique 
emphasized by Chen. For a distance d on N, we denote by 'Pd(N) the set of proba- 
bility measures £ on N such that X^eN d(x, x )£(x) < oo for some (or equivalently for 
all) Xq G N. We recall that the Wasserstein distance between two probability measures 
Pi, P>2 £ "Pd(N) is defined by 



W d (/!i,/! 2 )= inf / / d(x,y)>y(dx,dy), (3.1) 

76Marg(/ii,/i 2 ) JN JN 

where Marg(^i, P2) is the set of probability measures on N 2 such that the marginal 
distributions are pi and p2, respectively. The Kantorovich- Rubinstein duality [Vj th. 5.10] 
gives 

yV d {pi,P2) = sup gd(p 1 -p 2 ), (3.2) 
seLi Pl (d) « 
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where Lip(d) is the set of Lipschitz function g with respect to the distance d, i.e. 



\9\\up(d) 



sup 



\a(x) -g(y)\ 

d(x,y) 



< oo, 



and Lip 1 (d) consists of 1-Lipschitz functions. We assume that the kernel Ff G "Pd(N) for 
every ieN and t > so that the semigroup is well-defined on Lip(d). The Wasserstein 
curvature of (X t )t>o with respect to a given distance d is the optimal (largest) constant 
a in the following contraction inequality: 

t > 0. (3.3) 



I -ft 1 1 Lip(d)— >Lip(d) < e ° > 



Here ||Pt||Lip(d)->Lip(d) denotes the supremum of ||Pt/||Lip(d) when / runs over Lip 1 (d). It 
is actually equivalent to the property that 

W d {Ff i Ff)<e-' t d(x ) y) i x,y eN, t > 0. 

If the optimal constant is positive, then the process is positive recurrent and the semigroup 
converges exponentially fast in Wasserstein distance Wd to the stationary distribution p 
[Che2l th. 5.23]. 

Let p G T+ be an increasing function and define u G J 7 + as« x := p(x + 1) — p(x). The 
metric under consideration in the forthcoming analysis is 

d u (x,y) = \p(x) - p(y)\. 

Hence u remains for the distance between two consecutive points. In particular the space 
of functions / for which the intertwining relation of theorem 12.11 is available is actually 
Lip(du). Then it is shown in |Chel|. [J] by coupling arguments that the Wasserstein 
curvature a u with respect to the distance d u is given by the Chen exponent, i.e. 

a u = vat v x+ x - v x V \ x - K+\ • 

xeN u x u x 

The following corollary of theorem 12.11 allows to recover this result via an intertwining 
relation. 

Corollary 3.1 (Contraction and curvature). Assume that the potential V u is lower bounded. 
Then with the notations of theorem \2.1\ for any t > 0, 



■t||Lip(d u )-4-Lip(du) — 1 1 1 1 Lip(rf^) — SUp E x 



exp 



V U (X U S ) ds 



(3.4) 



In particular, the contraction inequality (13. 3p is satisfied with the optimal constant 

au = in f, v u{y)- ( 3 - 5 ) 

Proof. Let / G Lip 1 ((i u ) be a 1-Lipschitz function with respect to the distance d u . For 
any y, z G N such that y < z (without loss of generality), we have by the intertwining 
identity (12.21) of theorem 12.11 and Jensen's inequality, 

2-1 



\P t f(z)-P t f(y)\<J2u x \d u P t f(x)\ 



x=y 

2-1 



< J2 u * E * 



x=y 



\d u f{X U;t )\ exp 



V u (X UiS ) ds 



< d u (z,y) sup E x 

x&N 



exp (- J V u {X u>s )ds 
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so that dividing by d u (z,y) and taking suprema entail the inequality: 



I Pt 1 1 Li P (d u )->Lip(d u ) < sup E x 



exp 



V U (X U>S ) ds 



Finally, since by remark 12.41 the semigroup (Pt)t>o propagates monotonicity, the right- 
hand-side of the latter inequality is nothing but H-PfpHLip^), showing that the supremum 
over Lip 1 ((i u ) is attained for the function p. The proof of equation (I3.4|) is achieved. 

To establish (I3.5p . note that it suffices to get part < since the other inequality follows 
from f)3.4p . Applying (12.21) to the function p which is trivially in Lip 1 (c/ U ), we have for all 
x G N, 

1 



t 



log E x 



exp 



V u (X UiS ) ds 



t > 0, 

available for all iGN. 

□ 



and taking the limit as t — > entails the inequality o u < V u (x) 
The proof of (13.51) is now complete. 

Remark 3.2 (Pointwise gradient estimates for the Poisson equation). The argument 
used in the proof of corollary IJ.il allows also to obtain pointwise gradient estimates for 
the solution of the Poisson equation at the heart of Chen-Stein methods |B-H-Jl IBr-X} 
IB-Xt ISch] . More precisely, let us assume that d u is such that p G ^(p). For any 
centered function f G Lip 1 (c? u ), let us consider the Poisson equation —Cg = f , where 
the unknown is g. Then under the assumption a u > 0, there exists a unique centered 
solution gf G Lip(d u ) to this equation given by the formula gf = / °° P t f dt. We have for 
any i6N the following estimate (compare with |L-M| th. 2.1]/- 



sup 

/6Lip 1 (d u ) 



I % 0*0 1 



sup u 

/eLip^du) 



= Ux 
< —■ 

On, 



\d u P t f{x)\dt 
d u P t p(x) dt 



E, 



exp 



V U {X U S ) ds 



dt 



3.2. Functional inequalities. Theorems 12.11 and 12.71 allow to establish a whole family 
of discrete functional inequalities. We define the bilinear symmetric form T on J 7 by 



T(f,g):=-(C(fg)-fCg 



gCf) = -(Xdfdg + ud*fd*c 



Under the positive recurrence assumption, the associated Dirichlet form acting on its 
domain X>(£ M ) X T>(£,j) is given by 



£/i(/> 9) ■= \ I g)dp= [ \df dg dp 
2 Jn Jn 



where the second equality comes from the reversibility of the process. Here the domain 
T>{£^) corresponds to the subspace of functions / G L 2 (p) such that Sfj,{f,f) is finite. 
The stationary distribution p is said to satisfy the Poincare inequality with constant c if 
for any function / G T>{E^), 

cVar M (/) <£„(/,/), (3.6) 
where Var jU (/) := p(f 2 ) — p(f) 2 and p(f) := / N / dp. The optimal (largest) constant 
cp is the spectral gap of C, i.e. the first non-trivial eigenvalue of the operator —£. The 
constant cp governs the L 2 (p) exponential decay to the equilibrium of the semigroup: for 
all / G L 2 [p) and t > 0, 

ll^/-M/)IU 2M <e- cpt ||/-M/)IU^). 
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Several years ago, Chen used a coupling method which provides the following formula for 
the spectral gap: 

cp = sup o u 

where a u is the Wasserstein curvature of section 13.11 or, in other words, the Chen expo- 
nent. It corresponds to theorem 1.1 in |Chelj . equation (1.4). The following corollary of 
theorem 12.11 allows to recover the > part of Chen's formula. 

Corollary 3.3 (Spectral gap and Wasserstein curvatures). Assume that there exists some 
function u G J-+ such that the associated Wasserstein curvature a u is positive. Then the 
Poincare inequality ( 13.61) holds with constant sup M6 j- + a u , or in other words 

cp > sup a u . 

Proof. Since there exists some function u G J 7 + such that the Wasserstein curvature a u is 
positive, the process is positive recurrent. By proposition 6.59 in |Che2] . the subspace of 
T>{£^) consisting of functions with finite support is a core of the Dirichlet form and thus 
we can assume without loss of generality that / has finite support. We have 

r poo A 

Var M (/) = -/ / -(P t f) 2 dtdfi 
Jn Jo at 

= -2 / / P t fCP t fdtdfi 
Jn Jo 

roo p 

= 2 / Xu 2 (d u P t f) 2 dpdt 
Jo Jn 

poo p 

<2 / e- 2 ^ / Xu 2 P u , t (d u f) 2 d^dt, 
Jo Jn 

where in the last line we used theorem 12.11 with the convex function (p(x) = x 2 . Now the 
measure \u 2 fi is invariant for the semigroup (P u ,t)t>o, so that we have 

roo p 

Var M (/) < 2 / e- 2 ^ / \u 2 (d u f) 2 d^idt 
Jo Jn 

= — [ X(df) 2 df, 
a u Jn 

= —S^(f,f), 

where in the second line we used o~ u > 0. The proof of the Poincare inequality is complete. 

□ 

Remark 3.4 (M/M/oo and M/M/l). The spectral gap of the M/M/oo and M/M/l 
processes is well-known |Chel] . Corollary allows to recover it easily. Indeed, in the 
M/M/oo case, the value Cp = v can be obtained as follows: choose the constant weight 
u = 1 to get cp > v, and notice that the equality holds for affine functions. For a positive 
recurrent M/M/l process, i.e. X < v , we obtain cp > (v^A — \/v) 2 by choosing the weight 
u x '■= {ylX) x l 2 , whereas the equality asymptotically holds in f!3.6j) as k —> yvj\ for the 
functions G N. We conclude that cp = (s/X — \/v) 2 ■ 

Remark 3.5 (Alternative method for M/M/l). In the M/M/l case, let us recover the 
bound Cp > (v^A — y/v) 2 by using a different method. Letting p(x) := x for x G N and 
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g = f — /(0) for a given function f G V(£^), we have 

I 9 2 dfi = —^— I g 2 {-Cp)dp 
Jn v — A Jn 



< 



1 


V 


— 


A 




1 




V 




A 




A 




V 




A 




A 




V 




A 




A 




V 




A 



d(g 2 )dpdp 
J N (2gdf+(df) 2 ) dp 



g*dfiJ I \df)*dfi+ I (df) 2 dp 
V Jn Jn 



where in the last inequality we used Cauchy-Schwarz' inequality. Solving this polynomial 
of degree 2 entails the inequality 

A 



(df) 2 dp. 



Finally using the inequality Var At (/) < f N g 2 dp, we get the result. 

Remark 3.6 (Diffusion case). As mentioned in remark IK~8\ the argument above leading 
to the Poincare inequality might be extended to the positive recurrent diffusion case. In 
particular, under the same notation we obtain the following lower bound on the Poincare 
constant 

cp > sup inf V a (x), 

where the supremum is taken over all positive ff 2 function a on R. Note that up to the 
transformation a — > l/a, such a formula was already obtained by Chen and Wang in 
|C-W] through their theorem 3.1, equation (3.4), by using a coupling approach somewhat 
similar to that emphasized by Chen in the discrete case. 

Theorem 12.71 allows to derive functional inequalities more general than the Poincare 
inequality. Let X be an open interval of R and for a smooth convex function (p : X — > R 
such that ip" > and —l/(p" is convex on X, we define the (^-entropy of a sufficiently 
integrable function / : N — > X as 

Eat£(/) =/*(¥>(/))-¥> W))- 

Following |Chalj . we say that the stationary distribution p satisfies a ^-entropy inequality 
with constant c > if for any X- valued function / G T>{E^ such that <p'(f) G T>(£^), 

C Ent^(/)<^(/,^(/)). (3.7) 

See for instance [Cha2] for an investigation of the properties of (^-entropies. The (p- 
entropy inequality (13.71) is satisfied if and only if the following entropy dissipation of the 
semigroup holds: for any sufficiently integrable X- valued function / and every t > 0, 

Ent£(PJ)<e- ci Ent£(/). 

We have the following corollary of theorem 12. 7[ 

Corollary 3.7 (Entropic inequalities and Wasserstein curvature). If the birth rate A is 
non-increasing and the Wasserstein curvature a% (with the constant weight u = 1) is 
positive, then the if -entropy inequality (13.71) holds with constant o\. 
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Proof. As in the proof of corollary 13.31 the assertion o~% > entails the positive recurrence 
of the process. Moreover, we assume once again that the X-valued function / has finite 
support. By reversibility, we have 



Ent£(/) = / (<p(P f) - tpW))) dfx 
* Jn 

r roo d 

ip{Ptf) dtdfj, 



n Jo dt 

if'(Ptf) CPJdfidt 



JN 



JN 







\dP t fdip'(P t f)dfxdt 
XB^ (P t f,dP t f) dfidt, 



where is as in theorem 12.71 (the identity dg dip'(g) = B ip {g 1 dg) comes from g + dg = 
g(- + 1)). Using now theorem 12.71 together with the invariance of the measure A/z for the 
1-modification semigroup (Pi t t)t>o, we obtain 

Ent£(/) < / / e-^\P 1>t BV(f,df) d^dt 

* JO JN 

roo r 

= / / e- ait X B* (/, Of) dfxdt 

JO JN 

\B*(f,df) dfi 



G\ JN 

□ 

Remark 3.8 (Examples of entropic inequalities). The constant in the tp -entropy in- 
equality provided by corollary \3. 71 is not optimal in general (compare for instance with 
the Poincare inequality of corollary \3.3\ when (p(r) = r 2 with I = M.). The choice 
(p(r) = r logr with X = (0, oo) allows us to recover the modified \og-Sobolev inequality of 
[C-DP-P[ th. 3.1]: for any positive function f G T>[S^) such that log / G T>{£^), 

a.Ent^f) KS.if, log f). (3.8) 

Note that beyond this entropic inequality, it is proved in |C-DP-P] that the entropy is 
convex along the semigroup (a careful reading of the proof in |C-DP-P] suggests that it 
simply boils down to commutation and convexity of A transforms!). For the M/M/oo 
process, the estimate of corollary VS. 71 is sharp since o\ = v and the equality in (13. 8p holds 
as a — > oo for the function x G N e ax . Note that the M/M/l process and its invariant 
distribution, which is geometric, do not satisfy a modified log-Sobolev inequality. Another 
if -entropy inequality of interest is that obtained when considering the convex function 
<p(r) := r p , p G (1,2], with I = (0, oo): for any positive function f G T)[Ey) such that 

f" : < vu:,,). 

Kf p )-Kff<-£Mf p - 1 )- (3-9) 

Such an inequality has been studied in [B-T] in the case of Markov processes on a finite 
state space and also in |Cha2] for the M/M/oo queuing process. In particular, it can be 
seen as an interpolation between Poincare and modified log-Sobolev inequalities. 

Under the positive recurrence assumption, theorem 12. II implies also other type of func- 
tional inequalities such as discrete isoperimetry and transportation-information inequal- 
ities. Given a positive function u, we focus on the distance d u constructed in section I3~T1 
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where we assume moreover that p G X>(£ M ), i.e. Xu 2 is /i-integrable or, in other words, the 
^-modification process (X Ut t)t>o is positive recurrent. The invariant measure /i is said to 
satisfy a weighted isoperimetric inequality with weight u and constant h u > if for any 
absolutely continuous probability measure n with density / G T>{£y) with respect to /i, 

h u W du (7T,fx) < [ \u\df\dfi, (3.10) 

Jn 

where the Wasserstein distance Wd u is defined in (13.1 j) with respect to the distance d u . 
The terminology of isoperimetry is employed here because it is a generalization of the 
classical isoperimetry, which states that the centered L 1 -norm is dominated by an energy 
of L 1 -type. Indeed, if the weight u is identically 1, then the distance d\ between two 
different points is at least 1, so that 03. 101) entails 

2h 1 f \f-l\d{i = h 1 W d (n, f j)<h 1 W dl (n, f Ji l )< I X\df\dfi, 
Jn Jn 

where d is the trivial distance or 1. Note that the Z^-energy emphasized above differs 
from the discrete version of the diffusion case, since our discrete gradient does not derive 
from T. 

On the other hand, let us introduce the transportation-information inequalities em- 
phasized in [G-L-W-Y] . Let a be a continuous positive and increasing function on [0, oo) 
vanishing at 0. The invariant measure \i satisfies a transportation- information inequality 
with deviation function a if for any absolutely continuous probability measure n with 
density / with respect to /i, we have 

a(nU7r,^))<Z(7r,/i), (3.11) 

where the so-called Fisher-Donsker-Varadhan information of 7r with respect to fi is defined 

as 

^ ' 1 oo otherwise. 

Note that Z(-, \i) is nothing but the rate function governing the large deviation principle 
in large time of the empirical measure L t := t -1 Jq5x s ds, where 5 X is the Dirac mass 
at point x. In other words, the Fisher-Donsker-Varadhan information rewrites as the 
variational identity |Che2| th. 8.8]: 

, v r CV , 

X(tt, h)= sup / — — dn. 

The interest of the transportation-information inequality resides in the equivalence with 
the following tail estimate of the empirical measure [G-L-W-Y] th. 2.4]: for any absolutely 
continuous probability measure 7r with density / G L 2 (fi) with respect to /i, and any 
g G Lip^du), 

F„(L t (g)-ii(g) > r) < \\f\\ L 2 M e~ a ^ , r > 0, t > 0. 
We have the following corollary of theorem 12.11 

Corollary 3.9 (Weighted isoperimetry and transportation-information inequality). With 
the notations of theorem \2.1\ assume that the process is positive recurrent and that the 
following quantity is well-defined: 



Hi, 



[ sup E x exp (- [ V U (X U S ) ds) 

JO rceN L V Jo / 



dt < oo. 



Then the weighted isoperimetric inequality (I3.10P is satisfied with constant h u = 1/k u . 
If moreover there exists two constants e > and 9 > 1 such that 

(l + e)X x u 2 x + (l + l/e)u x ul_ 1 <-a(X x (9-l)+u x (l/9-l)) + b, xeN, (3.12) 
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where a := a e £ > and b := b £ _g > are two other constants depending on both e and 9, 
then the transportation-information inequality (13.111) is satisfied with deviation function 



Jb 2 + 2a(r/ Ku y-b 

a(r) := sup . 

e>o,0>i 2a 

Remark 3.10 (The case of positive Wasserstein curvature). In particular if the Wasser- 
stein curvature o~ u with respect to the distance d u is positive, then the process is positive 
recurrent and we have 

o- u W du {TT,fi) < \u\df\dfi and a (W du (-K, /i)) < I(ir, fi), 
Jn 

with the deviation function 



/ft 2 + 2a{ra u f - b 

a(r) : = sup . 

£>o,6>>i 2a 



Proof. For every f,g£ T^{£^) we have, by reversibility, 

Cov^(f,g) := J^(g - J^gdf?J f dp 



jf (- | o °° CP t gdt) fdfx 

r i~ L ptgiCfdfi ) dt 

S ll (P t g,f)dt. (3.13) 



Now, for every probability measure 7r <C fi with dn = fdfi, f G X>(£ M ), we get, using 

dsn, 

Wd u (tt, fi) = sup Gov M (/, g) 

gSLip 1 (d n ) 

TOO 



POO 

sup / £p(P t g, f) dt 
g eu Pl (d u ) Jo 

poo p 

sup / Xudf d u P t g djidt 

g£Li Pl (d u ) JO JN 



oo P 

/ \u\df\ d u P t pdjjdt 

JN 

t 



< 



/ sup exp ( — / V u (X us )ds) dt / \u\df\dfi 
Jo xeN L \ io J \ 



'0 xeN 

where in the last inequality we used theorem 12. II This concludes the proof of the weighted 
isoperimetric inequality. 

Using now Cauchy-Schwarz inequality, reversibility and then (I3.12p with Vg(x) := 9 X , 

x e N, 



2 

f ) dfi 



< k u fj.) y^T ((1 + e)\u* + (1 + l/e^u 2 ^) f d/i 

< k u ^X(tt, fi) J (-a + b^j fdfi 
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from which the desired transportation-information inequality holds. □ 

Remark 3.11 (M/M/oo and M/M/l revisited). Corollary \3.tA exhibits optimal func- 
tional inequalities, at least in the M/M/oo case and its stationary distribution, the Pois- 
son measure of mean X/u. Choosing the weight u = 1, we obtain the optimal constant 
hi = v in the isoperimetric inequality. Indeed, corollary VJ.iA entails h\ > v, whereas the 
other inequality is obtained by choosing it a Poisson measure of different parameter. For 
the transportation-information inequality, we recover theorem 2.1 in [M-W-W] since the 
choice of a := #(1 + l/e)/(6 — 1) and b := A(l + e + (1 + l/e)9) allows us to obtain the 

deviation function a(r) := A(y / 1 + ur/X — l) 2 , r > 0. Note that it is optimal in view 
of example 4-5 in |G-G-WJ : for any absolutely continuous probability measure 7r with 
square-integrable density with respect to \i, 



1 (l /•* 

lim - logPj - / X s ds 
t \t Jo 




, , ~o-„ , , , -~ ■ • , = -X(Jl + ^--l) . r>0. 

t \t Jo v ) V V A 



For the M/M/l process, we have the following inequalities for the optimal isoperimetric 
constant h u , with u x = (u/X) x ^ 2 (a quantity that will appear again in section UT^I : 

(Vx - v^) 2 < hu < - \fx)V^. 

To get the second inequality, we choose the density f = {y/X){l — l{o}) an d the 1-Lipschitz 
test function g = p. In particular as the ratio X/u is small, we obtain h u w v. However, 
we ignore if such a process satisfies a transportation-information inequality. 

3.3. Hitting time of the origin by the M/M/l process. Recall that we consider 
the ergodic M/M/l process (A < v) for which the stationary distribution is geometric 
of parameter X/v. Since the process behaves as a random walk outside 0, the ergodic 
property relies essentially on its behavior at point 0. Using the notation of theorem 12.11 
the intertwining relation (12.21) applied with a positive function u entails the identity 



d u PJ(x) = 



d u f(X t ) exp (- f\{X u , s )ds^ 



where the potential is given for every x G N by 

T/ ( \ u x-l ^ , u x+ i 
V u (x) :=v ^l{^o} + A A. 

Following Robert [Rj, the process (X^) t > is the solution of the stochastic differential 
equation 

XI = y and dXf = dN t W - l {x v_ >0} dN^ \ t > 0, (3.14) 

where (iV t (A) ) t > and (N t {u) ) t>o are two independent Poisson processes with parameter A 
and v, respectively. Since the process is assumed to be positive recurrent, the hitting 
time of 0, 

T y := inf{t > : X\ = 0} 
is finite almost surely. We have the following corollary of theorem 12.11 

Corollary 3.12 (Hitting time of the origin for the ergodic M/M/l process). Given 
x E N, consider a positive recurrent M/M/l process (Xf +1 ) t >o starting at point x + 1, 
and denote (X* t ) t > its u-modification process starting at point x, where 

«.:= (£)'>!. 
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Then we have the following tail estimate: for any t > 0, 

P(T* +1 >t)= u x e-< (^^) 2 E -^L-j exp (-VX^ J*l {0 y(X^) ds 



< u T e 



-^(^A-^/I7) ^ 



Proof. Let us use a coupling argument. Let (Xf ) 4 >o be a copy of (Xf +1 ) t >o, starting at 
point x. We assume that it constructed with respect to the same driving Poisson pro- 
cesses [Nf)t>o and (N^) t > as the process (Xf +1 ) t > . Hence the stochastic differential 
equation (I3.14|) satisfied by the two coupling processes entails that the difference between 
(X? +1 ) t > and (Xf) t > remains constant, equal to 1, until time Tq +1 , the first hitting 
time of the origin by (Xf +1 ) t > . After time Tq , the processes are identically the same, 
so that the following identity holds: 



x+l 



X^ + lrj,x+l 



{T* +i >t}' 



t > 0. 



Since the original process is assumed to be positive recurrent, the coupling is successful, 
i.e. the coupling time is finite almost surely. Therefore we have for any function / £ 
Lip(di), where d\ is the distance di(x,y) = \x — y\, 



/(*, 



dPJ(x) = PJ(x + 1) - PJ(x) = E 
so that if we denote the function p(x) = x, we obtain 

P(T* +1 >t) = dP t p(x) = u x d u P tP (x) 
Using now (\2.2\i with the function it, we get 

1 



E 



Bf(Xf) lr T »+l 



{Tg +l >t} 



F(T* +L > t) = u x E 



exp 



Vu(X*)d8 



where V u := - ^/vf + y/Xu 1 



{0}- 



□ 



Remark 3.13 (Sharpness). Using a completely different approach, Van Doom established 
in [VDj . through his theorem 4-2 together with his example 5, the following asymptotics 

lim - log P(7T +1 >t) = - v^7) 2 , x e N. 

Hence one deduces that the exponential decay in the result of corollary is sharp. On 
the other hand, proposition 5.4 in [R) states that Tq' +1 has exponential moment bounded 
as follows: 



E 



< 



A 



so that Chebyshev's inequality yields a tail estimate somewhat similar to ours - although 
with a worst constant depending on the initial point x + l. 

Remark 3.14 (Other approach). The proof of corollary suggests also a martingale 
approach. First, note that we have the identity 
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which entails as in the previous proof and since u > 1, the following computations: 
P(T* +1 >t) = dP tP {x) 

exp \ - J ul {0} (X x )ds 
' Cu 



= E 
< E 



< ^ e 



^(v'A-v^) 2 



since t/ie process (M^) t >o given by 

M? := exp 
z's a supermartingale. Indeed, denoting 



\ —(x:)ds), t>o, 

Jo u J 



Z u t :- 



exp 



* Cu 



u 



(x:) ds 



we have by Ito's formula: 

dM? = Z^du{X x ) + u[X x t ) dZ? 



Cu 



= Z^(dM t + Cu(Xf)dt) - u(X x ) — (XT) Z u t dt 

u 

= Z u t dM u 

where (M t )t>o is a local martingale. Therefore, the process (M") t >o is a positive local 
martingale and thus a supermartingale. 



3.4. Convex domination of birth-death processes. Let (Xf) t > be the M/M/oo 

process starting from x G N. The Mehler-type formula (II. 2p states that the random 
variable X x has the same distribution as the independent sum of the variable X t °, which 
follows the Poisson distribution of parameter (A/V)(l — e~ ut ), and a binomial random 

variable B[ x ^ of parameters (x, e~ ut ) 



By convention, is assumed to be 0. Hence we 



have for any non-negative function / and any iGN, 

E[f(X?))=E[f(X? + BM)), 



t > 0. 



(3.15) 



Such an identity can be provided by using the commutation relation (11.41) . Indeed we 
have 



E 



1 - e-*) E [f(X x )] + e"* E[f(X x + 1)] 



so that a recursive argument on the initial state provides the required result. An inter- 
esting consequence of (13.151) appears in terms of concentration properties. For instance 
a straightforward computation entails that for any 9 > 0, we get the following inequality 
on the Laplace transforms 



E 



< E 



ON? 



where N x is a Poisson random variable with the same mean as X x . Therefore, using 
the exponential Chebyshev inequality entails an upper bound on the tail of the centered 
random variable Xf — E[Xf\, which is sharp as t — > oo (recall that the stationary distri- 
bution is Poisson with parameter X/v). 

Actually, one may ask if for a more general birth-death process, the intertwining relation 
of type (12.21) may imply a relation similar to (I3.15p . This leads to the notion of stochastic 
ordering. 
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Following the presentation enlighten by Stoyan in [S], let us start with the classical 
notion of stochastic ordering for integer-valued random variables. We say that X is 
stochastically smaller than Y, and we note X <d Y, if for any function / G J~d, 

E[f(X)]<E[f(Y)}. 

Such a relation, as the convex domination introduced below, is a partial ordering on the 
set of distribution functions. The interesting feature of this stochastic ordering resides 
in its characterization in terms of coupling: we have X <d Y if and only if there exist 
random variables X\ and Y\, both defined on the same probability space and with the 
same distribution as X and Y respectively, such that P(ATi < X 2 ) = 1. Moreover, it 
is equivalent to the following comparison between tails: we have X <d Y if and only 
¥(X > x) < P(V > x) for any x G R. In other words, the random variable X takes small 
values with a higher probability than Y does. 

Another stochastic ordering of interest is the convex ordering, or convex domination. 
Denote T c the subset of consisting of non-negative non- decreasing convex functions, 
where in our discrete setting the convexity of a function / : N — > R is understood as 
d 2 f > 0. We say that X is convex dominated by Y, and we note X < c Y, if for any 
function / G T c , 

E[f(X)]<E[f(Y)]. 
It is known to be equivalent to the inequality 



E 



{X-x)+ <E (Y-x 



x G 



where a + := max{a,0}. Typically, one may deduce from the convex domination concen- 
tration properties like a comparison of moments or Laplace transforms as in the M/M/oo 
case above. Moreover, this refined ordering might appear for instance when using de- 
la- Vallee-Poussin's lemma about uniform integrability of a family of random variables. 
However, in contrast to the <d ordering, the authors ignore if there exists a genuine 
interpretation of the convex domination in terms of coupling. 

Coming back to our birth-death framework, we observe that if we want to use the 
intertwining relation (12.21) of theorem 12.11 in order to obtain stochastic domination, then 
a first difficulty arises. Indeed, another birth-death process appears in the right-hand- 
side of (12.21) . namely the u- modification of the original process. Therefore, let us provide 
first a lemma which allows us to compare two birth-death processes with respect to the 
<d ordering. Although the result below is somewhat obvious from the point of view of 
coupling, we give an alternative proof based on the interpolation method emphasized in 
the proof of theorem 12.1 [ See also [HI prop. 4.2.10]. 

Lemma 3.15 (Stochastic comparison of birth-death processes). Let (Xf) t >o and (X£)t>o 
be two birth-death processes both starting from x G N. Denoting respectively A, v and A, v 
the transition rates of the associated generators £ and C, we assume that they satisfy the 
following assumption: 

A < A and v > v 

Then for every t > 0, the random variable Xf is stochastically smaller than X* . In other 
words, we have Xf <d Xf . 

Proof. Let g G J~d and define the function s G [0,t] H- J(s) := P s P t - s g where (Pt)t>o an d 
(Pt)t>o ar e the semigroups of (Xf) t > and {Xf) t >o respectively. By differentiation, we 
have 

J'(s) = P s (£P t _ s g - CP t _ s g) = P s ((A - A) dP t „ s g + {u - v) d*P t _ s g 
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which is non-positive since the semigroup (Pt)t>o satisfies the propagation of monotonicity, 
cf. remark 12.41 Hence the function J is non-increasing and the desired result holds. □ 

Now we are able to state the following corollary of theorem 12. 11 which states a new con- 
vex domination involving decoupled random variables in the right-hand- side. However, 
despite some particular cases like the M/M/l case for which the convenient coupling 
appearing in the proof of corollary 13.121 allows us to extend the next result to the <^ 
ordering, we ignore if it can be done in full generality. 

Corollary 3.16 (Convex domination). Denote (Xf) t > a birth-death process starting at 
some point y G N. We assume that the birth rate A is non-increasing and that there exists 
K > such that 

d{y - A) > k. 

Then for any t > and any x G N, the random variable Xf +1 is convex dominated by the 
independent sum of Xf and a Bernoulli random variable Y t of parameter e~ Kt G (0,1]. 
In other words, we have 

X? +1 < c Xf + Y t . 

Proof. We have to show that for any function / G J- c , 

E[/pO] <E[/(.x? + y t )]- (3.i6) 

Using the intertwining relation (12.21) of theorem 12.11 we have since / is non-decreasing: 



E 



f(Xr i ) < E[/(Xf )] + e~ E df{X* 



<E[f(Xn}+e- Kt E[df(Xn) 

= (1 - e~ Kt ) E[/(Xf )] + e- Kt E[f(X? + 1)] 

= E\f(Xf + Y t )], 

where to obtain the second inequality we used lemma 15. 15l with the 1-modification process 
{Xf t ) t >o playing the role of (Xf) t >o since df is non-decreasing (recall that / G J- c ). □ 

Remark 3.17 (More on convex domination). By an easy recursive argument one obtains 
from the latter result the following convex domination: 

X?< c X? + Bt\ 

where B^ is a binomial random variable of parameters (x,e~ Kt ), independent from X® , 
as in the case of the M/M/oo queuing process. 
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